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Abstract. In this paper we study the localization transition induced by the disorder in random antiferro- 
magnetic spin-i chains. The results of numerical large scale computations are presented for the XX model 
using its free fermions representation. The scaling behavior of the spin stiffness is investigated for various 
disorder strengths. The disorder dependence of the localization length is studied and a comparison between 
numerical results and bosonization arguments is presented. A non trivial connection between localization 
effects and the crossover from the pure XX fixed point to the infinite randomness fixed point is pointed 
out. 
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1 Introduction 

Quantum spin chains exhibit a large number of interesting 
features because the quantum fluctuations are often rel- 
evant, especially at low temperature. The antiferromag- 
netic (AF) Heisenberg model in one dimension (ID) has 
been extensively studied since the discovery in 1931 of the 
Bethe Ansatz pQ for the spin S = h chain. In ID, the AF 
XXZ model defined by the Hamiltonian 



where eo is the GS energy per site. It measures the mag- 
netization transport along the ring and in the fermionic 
language this is called the charge stiffness, which is the 
Drudc weight of the conductivity. The gap-less phase is 
characterized by peculiar transport properties: in the ther- 
modynamic limit Shastry and Sutherland |S] showed that 
in the critical regime the spin stiffness of the XXZ chain 
follows: 



H 



XXZ 



Ps(A) = J 



Si +1 + h.c) + AS?S, 



(1) 



7r sin(/x) 

4/i(7T - fi) 



where A — cos(p), (4) 



with J > and A > 0, exhibits a gap-less excitation spec- 
trum for S = \ if A < 1, whereas a gap opens up in the 
spectrum when A > 1. In ID, the quantum fluctuations 
prevent the formation of true long-range order [2] but in 
the critical regime A < 1 the model [Eq.lJTJ] displays a 
gwasz-long-range order (QLRO) with power-law decaying 
spin-spin correlation functions in the ground state (GS). It 
is well known that the model [Eq.QJ], without quenched 
disorder, is integrable for conventional periodic boundary 
conditions pP as well as in the more general case of twisted 
boundary conditions (TBC) [3]. The latter are defined by: 
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where <p is the twist angle and is equivalent to a ring of in- 
teracting fermions threaded by a magnetic flux of strength 
— <j> |3]. The spin stiffness ps is defined by 



Ps 



L 



2 d 2 eo(0) 



10=0, 



(3) 



and it vanishes for A > 1. The phase transition which 
occurs at A = 1 can be viewed as a metal-insulator tran- 
sition [S] between a critical metallic phase with a finite ps 
and a gaped insulating regime where ps = 0, following a 
Mott mechanism. 

When the system is not homogeneous, the situation 
described above changes dramatically. For instance when 
only one coupling exchange is weaker than the others in an 
otherwise homogeneous ring, the stiffness has been found 
to scale to zero by numerical studies 6 , in perfect agree- 
ment with renormalization group arguments of Eggert and 
Affleck [7], and Kane and Fisher 

Moreover, for the case of a random spin-^ chain, Doty 
and Fisher [S] performed a bosonization study consider- 
ing several types of random perturbations added to JIJ. 
They concluded that in the AF critical regime the GS with 
QLRO is destroyed by any small amount of disorder and 
the phase transition associated is an Anderson localization 
transition |ll)j . reminiscent of the localization problem in 
ID disordered metals studied by Giamarchi and Schulz 
[TT] . A relevant length scale associated with the Anderson 
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transition is the localization length £* . 
More generally, the problem of transport in ID random 
media ^5] as well as localization effects and persistent cur- 
rents in disordered quantum rings have motivated a large 
number of theoretical studies in recent years [161117111811191 
. In the context of mesoscopic physics it turned out 
to be very interesting to study the transport properties for 
finite systems, where coherence effects are important |22l 
I23| . In particular the finite size (FS) dependence of the 
current, susceptibility and stiffness are important for a 
complete understanding of the experimental results. 

In the present paper we investigate the scaling behav- 
ior of the spin stiffness of the random spin-i chain. It 
is organized as follows. In Sec. 2, the numerical method, 
based on the free fermions formalism, is explained and 
notably the computation of the spin stiffness is described. 
Sec. 3 is devoted to the study of the localization transi- 
tion: Using some bosonization arguments as well as FS 
scaling analysis, an universal scaling of the stiffness to 
is expected and we compare it with numerical results. In 
Sec. 4, the disorder dependence of the localization length 
is studied and the bosonization predictions are demon- 
strated to be valid only for weak randomness. For strong 
disorder we propose a new quantity which gives a better 
description for the disorder dependence of £* . The relation 
to crossover effects observed recently for spin-spin correla- 
tion functions |12II13| is also worked. Sec. 5 contains some 
concluding remarks. 



The sign of the boundary term depends on the parity of 
the total number of fermions Af = 53;=i ^Yjj indeed 



M0) 



^{Clde^ + C\C L e l4> ). (8) 



Hence, when <fi — the resulting free fermions problem 
must have anti-periodic boundary conditions if the num- 
ber of fermions is even and periodic boundary conditions 
if N is odd [21] • In the non-random case, the solution of 
the problem via a Fourier transformation is trivial j^HI due 
to its translational invariance. In fc-space, the pure model 
is given by 



H 



xx 



-jJ2 c l c kCos(k). 



(9) 



Its GS is at half- filling (Af = -|, corresponding to the 
S* ot = sector) . The twist angle at the boundary produces 
a shift in the momentum space k — > k + (f> which can be 
uniformly distributed over all bond resulting in a local 
twist 8(f> = for each bond. Therefore the GS energy per 
site takes the following simple expression 



zq(L,4>) = -j ^cos(- 2?n ' 



L 



Jcos(f) 

L } L S in(|) (U) 



from which we can easily extract the spin stiffness 26 



2 Numerical method at the XX point 

We start with the ID random XX model on a ring closed 
with TBC. It is defined by 



L-l 



n™ dom (<t>) - E [y 5 r+i + M] + Mfl. (5) 



with the boundary term /i.l(</>) 



Jl(q+ 



(Sjftf-e-* + h.c). 



The couplings Ji are independent random numbers. 



2.1 Free fermions formulation 

For S = 4 the well known Jordan- Wigner mapping trans- 
forms spin operators into Fermi operators according to 



K^Cj. (6) 
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Nj = CjCj is the occupation number (0 or 1) at site 
j, given by Nj = 1/2 + 5|. Note that the Fermi anti- 
commutation relations are satisfied {C\,Cj} = S ld . The 
Hamiltonian J5J can then be written as 

L-l j 

n?* dom w = E [7(^+1 + c l+i c i)} + M# (7) 

i=l 



p s (L) = JOLsiny)- 1 ~ 7T- 1 + 0(L- 2 ). (11) 

When the system is inhomogeneous, the translational in- 
variance is broken and a solution in the reciprocal space 
is no longer possible. Fortunately the problem can be eas- 
ily diagonalized numerically, using standard linear algebra 
routines |27| . Indeed, with an unitary transformation the 
hamiltonian (7J can be expressed in a diagonal form [251 
I24U28U29) . For completeness we give here a brief descrip- 
tion of the method. First let us define a column vector <P 
of size L and its conjugate row vector •Z'T by 



(12) 



Hence, using this notation, we can re-write the Hamilto- 
nian JJJ) in terms of a symmetric L x L band matrix A{4>) 
as 



n 



XX 
random 



{</)) = ^ A(<f>) W , 



(13) 



with non-zero elements given by = and at the 

boundaries, A\,l — {—l) J ^^-e~ % ^. One can define the uni- 
tary transformation P that diagonalizes A. Then we get 
a new set of Fermi operators r\ q defined by 



^ = E^, ?4 = E p V^ 



(14) 



which yields the following diagonal form for the Hamilto- 
nian 



9=1 



(15) 
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where the e q (4>) are the eigenvalues of A((f>). At temper- 
ature T, the occupation number is given by the Fermi 
function (N q ) = (77]?^) = (1 + e e «w/ r ) _1 . Because of 
the particle-hole symmetry, the eigenvalues occur in pairs, 
equal in magnitude and opposite in sign. Hence, at T = 0, 
the GS energy is simply given by 



M=L/2 

e oO) = e 9^)' 

9=1 

where ei(<f») < e 2 ((f>) < ■■■ < e L (4>)- 



2.2 Numerical evaluation of the spin stiffness 
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figure 1. (a) Magnitude of the FS corrections of the spin stiff- 
ness \ps(L) — 7r _1 | for different choices of the twist angle 4> 
calculated for the pure XX model [Eq.©]. The long-dashed 
line is the exact result Lsin(^)" 1 — 7r _1 and the different sym- 
bols show the numerical results for different values of the twist, 
(b) Function Y(8<f>) = 2 — £^ry^ computed in double precision 
type. 



Numerical estimates for the spin stiffness can be ob- 
tained by approximating Eq.(PJ for finite L by 



PS 



, eoW-eo(0) 



(17) 



where 5<fi = <fi/L is the twist per site. Hence for a given 
system the calculation of ps requires to compute Eq. (|16fl 
twice: once for finite <fi and once for <fi = 0. Since the correc- 
tions are of order 1/L 2 an extrapolation L — > oo is, in prin- 
ciple, straightforward and yields the desired result. How- 
ever, the appropriate choice of <f> is somewhat delicate as 
we show in FigQ](a). Here the numerical results for the FS 
scaling of the spin stiffness of the pure chain are depicted, 
computed for various system sizes {L — 4, 8, 16, 2048) 
with three different values of the twist angle, and com- 
pared to the exactly known result given by Ea. l|ll|) . The 



discrepancy between the numerical data and the exact re- 
sult, observed for Stp/n = 10~ 3 and Sip/n = 10~ 5 can 
be understood as follows. Using Eg. (fTTjjl one can rewrite 
Eq.(jni) as 



1 — COs(<5</>) T ,, . 7T. -i 



(18) 



The function Y(5(f) — 2 1 ^jf^ , which is exactly equal to 
one when 5(f> = 0, is expected to decrease slowly when 8<j) 
increases. However, the numerical calculation of Y(S(f>) is 
limited by the machine precision and therefore we observe 
in Fig^(b) that even in double precision type, for < 
10 -4 undesirable oscillations appear. This puts a bound 
for the smallest value of S(f> that is meaningful for our 
numerical procedure, as we demonstrate in FigQ (a) for 
S4>/tt = 1CT 5 . On the other hand, when S(f> > 1CT 4 the 
value of Y deviates significantly from one as shown in 
Figffl ( a ) f° r S(p/n = 10~ 3 . Therefore, for a numerical 
calculation in double precision, the numerical derivative 
Ea. (|17l) gives the most reliable results for 5<fi ~ 10~ 4 7r 
which is confirmed by the numerical data obtained in this 
case, shown in Fig^](a). 



3 Localization transition : Scaling from pure 
to infinite randomness behavior 

3.1 Bosonization predictions for weak disorder and 
scaling argument in the localized-random singlet phase 

The critical behavior of the XXZ model [Eq. with weak 
randomness in the couplings and/or in external magnetic 
fields has been studied by Doty and Fisher [5] using a 
bosonization approach. They found that for random per- 
turbations which preserve the XY symmetry, the critical 
properties belong to the universality class of Giamarchi- 
Schulz transition for ID bosons in a random potential |llj. 
Let us define the disorder parameter T> by 



V={Jif- (Ji)' 



(19) 



For weak initial randomness 2?o <ti 1, the renormalization 
of the disorder under a change of length scale I — In L is 

^ = (3-2*)*,, (20) 
where K is the Z\-dependent Luttinger liquid parameter 



K(A) 



Therefore, if K < 3/2 (i.e. -\<A<1) 



2(tt-m) ' 

the disorder is a relevant perturbation and the line of pure 
fixed points is unstable under any amount of randomness. 
Under renormalization the system runs into an infinite 
randomness fixed point (IRFP) |9"ll3()| . Using a real space 
decimation procedure |2Q|, Fisher reached the same con- 
clusion and demonstrated analytically the existence of an 
attractive IRFP. Strictly speaking, at the IRFP the sys- 
tem is in the so-called random singlet phase (RSP) or 
in the fermionic language, the fermions are localized and 
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their transport properties are the ones of an insulator. 
For instance, the Drude weight is expected to be in the 
thermodynamic limit L — > oo. The renormalization flow 
is controlled by a disorder dependent length scale which 
emerges from Ea. (|20[) . the localization length: 

C(V)~V~^?. (21) 

In the thermodynamic limit the spin stiffness is finite in 
the QLRO phase (see Eq.QJ) and its FS scaling behavior 
is well known On the other hand, when V > we 
have ps(L,T>) — > and expect a scaling of the form 

ps (L,V) = g{-^), (22) 

with g a universal function. Defining x — L /£*(!>), one 
can consider 3 different regimes : (i) For x <C 1, i.e. on 
small length scales, the systems appears to be delocal- 
ized with g ~ ir~ 1 . (ii) For x 3> 1, i.e. on large length 
scales, the system shows the asymptotic behavior of the 
IRFP and g — ► 0. (iii) In the intermediate region x ~ 1, 
a crossover between the pure repulsive fixed point and 
the attractive IRFP occurs. Utilizing standard FS scal- 
ing arguments |31|. one can predict the behavior of g(x) 
in the asymptotic regime of the IRFP: ps has dimension 
of inverse (length^ -2 X £ T ), where £ r is the correlation 
length in the imaginary time direction 31 . In our case 
£ T ~ exp(^ 1 / 2 ), which is one manifestation of the crit- 
ical behavior at the IRFP (i.e. the dynamical exponent 
formally being z = oo), and £ = L for a finite system at 
criticality. Therefore we expect ps to scale as 

In PS {L) ~ -VI. (23) 

Combining this with Ea. l|22|l . we expect g(x) to behave 
as a constant ~ tt~ 1 in the delocalized regime (i) and to 
vanish as 

lng(x) Vx (24) 

in the localized regime (ii). 

3.2 Numerical results 

Following the method explained in Sec. 2., we study the 
spin- i XX model [Eq.©] with random bonds Ji distributed 
according to the flat distribution 

V(J) = {^' li ' e[1 - W > 1 + W] , (25) 
y ' \ otherwise. v ; 

which implies that the disorder strength is V = -| W 2 . Due 
to the strong sample-to-sample fluctuations that occur in 
many disordered quantum systems at low or zero temper- 
atures we have to perform a disorder average over a suf- 
ficiently large number of samples. In our calculations the 
latter ranges from N s = 10 3 for the biggest size up to 10 5 
for the smaller ones such that the error bars are well con- 
trolled, as we checked carefully. The system sizes vary from 
Lmin = 8 to L max = 2048 and we considered a large range 



•— • L=8 
■ ■■■ L=16 
♦--» L=32 




w 

figure 2. Disorder averaged value of the spin stiffness ps vs 
the disorder strength W ~ y/V for different system sizes, as 
indicated on the plot. Averaging has been done over N s = 10 3 
samples for the biggest sizes and up to 10 s for the smallest 
ones such that the error bars are well controlled, as we can 
observe. The expected behavior in the thermodynamic limit is 
represented by the black stars. 



of disorder strengths between W m i n = 0.025 and W max — 
1. The spin stiffness ps was evaluated using Eq. lfTTjl with a 
twist angle (j> = LxStfi = Lx 7r/10 4 and was then averaged 
over N s independent samples: J§ = ^ J] {samples} ps- In 

FigI2ps(IT) is shown for different system sizes and we see 
clearly that it approaches zero for increasing L. In order to 
validate the FS scaling form [Ea. lf^3*|l ]. we studied the dis- 
tribution of Imps- For W = 0.5, FigOH(a) shows such a dis- 
tribution for system sizes ranging from 8 to 512 sites with 
N s = 10 4 samples. As expected for a system described 
by an IRFP the distribution gets broader with increasing 
system size, which confirms that the dynamical exponent 
is formally infinite z = oo. Following Ea. 1)23(1 . the distri- 
bution P(^£S-) is plotted in Fig|3] (b) and as expected, 
the data of Fig|3] (a) collapse in a universal function. 

When the disorder is weaker, we expect strong FS ef- 
fects and a disorder-dependent length scale might control 
a crossover between the pure repulsive XX fixed point 
and the attractiv e IRFP . Such a behavior is illustrated 
in Fig^a) since ps(L) has been calculated for various 
values of the disorder W . Typically, when W > 0.3 we 
can observe the asymptotic behavior In ps(L) ~ —L 1 / 2 as 
soon as L ~ 100 but when W < 0.1 the pure behavior 
~p~s — 7T _1 remains robust up to very large L and even for 
L = 2048 the IRFP asymptotic regime is not yet reached. 

In order to characterize this crossover behavior, we 
studied the scaling function defined by Eq. 1221) and a cor- 
responding scaling plot of — (In ^(L/^*)) -1 is shown in 
FigH». For W = 0.225 we have chosen £* = 100 such 
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figure 3. Distribution of the spin stiffness obtained with N s = 
10 4 samples at W = 0.5. The system sizes are indicated on the 
plot, (a) The distribution of ln(ps) is broadening with system 
size, (b) Scaling plot of the data shown in Fig[3](a), assuming 
that the logarithm of the stiffness varies as the square root of 
the system size. 



that the crossover region is centered around x — L/t;* ~ 1 
and the other estimates, indicated on the plot, have been 
adjusted carefully in order to obtain the best data col- 
lapse. The 3 regimes mentioned above (see Sec. 3.1) are 
clearly visible : The pure regime (i) for which the stiff- 
ness takes values close to 7r _1 is observed if x <C 1. When 
i> 1 the infinite randomness regime (ii) is relevant : the 
universality of the IRFP is recovered and g(x) is in per- 
fect agreement with Ea.l|24|l. The intermediate crossover 
regime (hi) is visible for x ~ 1. 



4 The localization length as a crossover 
length scale 

Finally, we study the disorder dependence of the localiza- 
tion length £*. Using the values extracted from the data 
collapse shown in FigBJb), £*(T>) is shown in FigOJa) for 
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figure 4. (a) Inverse logarithm of the disorder averaged spin 
stiffness plotted for several box sizes W specified on the plot. 
The error bars are smaller than symbol sizes. The full line 
stands for the pure case and the dotted one shows the expected 
IRFP behavior [Eo.ffity]. (b) Scaling plot according to Eq.J23 
of the data shown in Fig^Ia) with indicated on the plot for 
each W . Pure and IRFP behavior are indicated respectively by 
full and dotted lines. 



several values of the disorder strength. The numerical re- 
sults are compared with the predicted power-law behav- 
ior Eu.pT)l which is at the XX point given by £*(£>) ~ 
T> , The agreement between the numerical results and 
the bosonization prediction is very good for weak disor- 
der, but for V > 0.1 the data deviate from a power-law. 
In order to extract a functional form for £* also in this 
range of disorder we look at its behavior as a function of 
the variance (5 of the random variable In Jj : 



5 = 



\j{\^Ji) 2 - (in J, 



which is related to W via 



1 - W 2 
4W 2 



In 



1 + VF 
1 - W 



(26) 



(27) 
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figure 5. Disorder dependence of the localization length £* of 
the random XX chain. Numerical results are shown with open 
circles and full lines represent power-laws as indicated on the 
plot, (a) As a function of the disorder parameter T> and (b) as 
a function of the disorder parameter 5. 



As we can observe in FigGJb), the parameter S is very 
useful to describe the disorder dependence of £* for any 
strength of randomness, indeed the power-law £* (S) ~ <5~ 2 
works perfectly for the whole range of randomness consid- 
ered here. Hence we assume that Ea. (l21ll has to be re- 
placed, for strong disorder, by 

C(8)~S-*, (28) 

and since for weak disorder S ~ we expect <S> — 3 J* 2 k ■ 
Actually, a similar conclusion was drawn in |12II13| . 
where the crossover effects visible in the spin-spin corre- 
lation function of random AF spin chains were studied. 
Indeed the correlation functions of the weakly disordered 
spin-i XXZ chain display a strong crossover behavior con- 
trolled by a disorder-dependent crossover length scale £ 
which behaves as J- 1 - 8 * - 2 ^2]. In analogy to what we 
did with the stiffness above, we can extract the crossover 
length scale £ using the scaling function 

~c{x) = C avg (L)/Co(L), with x = j, (29) 

where Cq(L) and C avg {L) are spin-spin correlation func- 
tions calculated at mid-chain respectively for the pure 
and random models. At the XX point, when W = 0, 
C (L) oc L- 1 ' 2 and at the IRFP C avg (L) oc L~ 2 . The 
crossover between these two distinct behaviors is shown 
in FigUJa) where c(x) presents a universal form, following 
c(x) ^constant for x -C 1 and c{x) ~ £~ 3 / 2 for x 3> 1. 
We see that the characteristic length scale £ beyond which 
the asymptotic IRFP behavior sets in in the correlation 
function scales with disorder strength in very much the 
same way as the localization length £*. 

5 Conclusion 

In this paper we have studied the scaling behavior of 
the stiffness of the random antiferromagnetic spin-i XX 




figure 6. (a) Scaling plot according to Eg. 1291 of mid-chain 
xx correlation function data obtained in I I2ll.'-j| for 5 differ- 
ent values of W indicated on the figure as well as the £ used 
for the data collapse. The line with open circles shows the 
pure behavior and the full line shows the RSP behavior at the 
IRFP. The crossover length scale £ is plotted vs T> and fitted 
by x> -1 ' 048±0 - 1 only for weak disorder in (b) whereas in (c) £ (5) 
displays a better agreement with a power law ~ g-is±o.2^ 



chain numerically via exact diagonalization calculations 
utilizing the fact that the system can be mapped on a 
free fermion model. The latter allowed us to study rather 
large system sizes by which we were able to analyze thor- 
oughly the crossover effects observable for weak disorder. 
Our results clearly show that the asymptotic behavior of 
the model under consideration is governed by an infinite 
randomness fixed point for all disorder strengths, includ- 
ing the weakest, as predicted by D. Fisher PU]. We could 
observe one of the characteristics of the IRFP, namely a 
formally infinite value for the dynamical exponent, from 
the finite size scaling behavior of the probability distribu- 
tion of the stiffness, where hips/L 1 / 2 occurs as a scaling 
variable indicating that the stiffness scales exponentially 
with the the square root of the system size. 

Moreover we showed that the finite size scaling form 
of the average value of the stiffness is governed by a char- 
acteristic length scale that depends on the strength of the 
disorder. The length scale can be identified as a localiza- 
tion length with regard to transport properties but also 
as a crossover length scale below which the system be- 
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haves essentially like a pure (disorder free) chain and the 
stiffness is constant and beyond which the asymptotic be- 
havior characteristic for an infinite randomness fixed point 
becomes visible and the stiffness scales to zero with a char- 
acteristic power of the system size. We found that this 
length scale diverges like 1/8 2 with decreasing variance 
S of the disorder, which agrees well with an analytical 
prediction using bosonization techniques. This behavior 
agrees also well with the scaling behavior of the crossover 
length for the spin-spin correlation function, which indi- 
cates that there is indeed a single disorder strength de- 
pendent length scale governing the crossover as well as 
the localization phenomena in this system. 



References 

f. H. Bethe , Z. Physik 71, 205 (1931). 

2. N. D. Mermin and H. Wagner, Phys. Rev. Lett. 17, 1133 
(1966). 

3. C. J. Hamer, G. R. W. Quispel, and M. T. Batchelor, 
J. Phys. A 20, 5677 (1987); F. Woynarovich and H.-P. Eckle, 
J. Phys. A 20, L97 (1987); F. C. Alcaraz, M. N. Barber, and 
M. T. Batchelor, Ann. Phys. 182, 280 (1988). 

4. N. Byers and C. N. Yang, Phys. Rev. Lett. 7 46 (1961). 

5. B. S. Shastry, and B. Sutherland, Phys. Rev. Lett. 65, 243 
(1990); B. Sutherland and B. S. Shastry, Phys. Rev. Lett. 
65, 1833 (1990). 

6. T.M.R. Byrnes, R.J. Bursill, H.-P. Eckle, C.J. Hamer, and 
A. W. Sandvik, Phys. Rev. B 66, f953f3 (2002). 

7. S. Eggert and I. Affleck, Phys. Rev. B 46, f0866 (f992); S. 
Eggert and L Affleck, Phys. Rev. Lett. 75, 934 (f995). 

8. C. L. Kane and M. P. A. Fisher, Phys. Rev. Lett. 68, f220 
(f992); 46, 15 233 (1992). 

9. C. A. Doty and D. S. Fisher, Phys. Rev. B 45, 2164 (f992). 
fO. E. Abrahams, P. W. Anderson, D. C. Licciardello and T. 

V. Ramakrishnan, Phys. Rev. Lett. 42, 673 (f979). 
ff. T. Giamarchi and H. J. Schulz, Phys. Rev. B 37, 325 
(1988). 

12. N. Laflorencie and H. Rieger, Phys. Rev. Lett. 91, 22970f 
(2003). 

f 3. N. Laflorencie, H. Rieger, A. W. Sandvik, and P. Henelius, 

preprint cond-mat/03f 2572 
f4. K. Hamacher, J. Stolze, and W. Wenzel, Phys. Rev. Lett. 

89, 127202 (2002). 

15. O. Montrunich, K. Damle, and D. Huse, Phys. Rev. B 63, 
134424 (2001). 

16. G. Bouzerar, D. Poilblanc, and G. Montambaux, Phys. 
Rev. B 49, 8258 (f994). 

f7. K. J. Runge and G. T. Zimanyi, Phys. Rev. B 49, f52f2 
(f994). 

18. T. Giamarchi and B. S. Shastry, Phys. Rev. B 51, 10915 
(1995). 

19. P. Schmitteckert et al., Phys. Rev. Lett. 80, 560 (1998). 

20. L. Urba and A. Rosengren, Phys. Rev. B 67, 104406 
(2003). 

21. F. Schiitz, M. Kollar, and P. Kopietz, Phys. Rev. Lett. 91, 
Of 7205 (2003). 

22. L Affleck and P. Simon, Phys. Rev. Lett. 86, 2854 (2001); 
P. Simon and L Affleck , Phys. Rev. Lett. 89, 206602 (2002). 

23. W. G. van der Wiel et al., Science 289, 2f05 (2000); J. 
Nygard, D. H. Cobden, and P. E. Lindelof, Nature (London) 
408, 342 (2000). 

24. A. P. Young and H. Rieger, Phys. Rev. B 53, 8486 (f996). 

25. E. Lieb, T. Schulz, and D. Mattis, Ann. Phys. (NY) 16, 
407 (1961). 

26. N. Laflorencie, S. Capponi and E. S. Sorensen, Eur. Phys. 
J. B 24, 77 (200f). 

27. The linear algebra package LAPACK has been used here. 

28. P. Henelius and S. M. Girvin, Phys. Rev. B 57, fl457 
(1998). 

29. F. Igloi, R. Juhasz, and H. Rieger, Phys. Rev. B 61, 11552 
(2000). 

30. D. S. Fisher, Phys. Rev. B 50, 3799 (1994). 

31. M. Wallin, E. S. S0rensen, S. M. Girvin, and A. P. Young, 
Phys. Rev. B 49, 12115 (1994). 

32. Strictly speaking there is a logaritmic correction but it is 
subdominant. 



